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We introduce the idea of actually cooling quantum systems by means of incoherent thermal light, hence giving 
rise to a counter-intuitive mechanism of "cooling by heating". In this effect, the mere incoherent occupation of 
a quantum mechanical mode serves as a trigger to enhance the coupling between other modes. This notion of 
effectively rendering states more coherent by driving with incoherent thermal quantum noise is applied here to 
the opto-mechanical setting, where this effect occurs most naturally. We discuss two ways of describing this 
situation, one of them making use of stochastic sampling of Gaussian quantum states with respect to stationary 
classical stochastic processes. The potential of experimentally demonstrating this counter-intuitive effect in 
opto-mechanical systems with present technology is sketched. 



Cooling in quantum physics is usually achieved in just the 
same way as it occurs in classical physics or in common ev- 
eryday situations: One brings a given system into contact with 
a colder bath. Coherent driving of quantum systems can effec- 
tively achieve the same aim, most prominently in instances of 
laser cooling of ions or in its opto-mechanical variant, cooling 
mechanical degrees of freedom using the radiation pressure of 
light. The coherence then serves a purpose of, in a way, ren- 
dering the state of the system "more quantum". In any case, in 
these situations, the interacting body should first and foremost 
be cold or coherent. 

In this work, we introduce a paradigm in which thermal 
hot states of light can be used to significantly cool down a 
quantum system. To be specific, we will focus on an opto- 
mechanical |Q}{6) implementation of this idea: This type of 
system seems to be an ideal candidate to demonstrate this ef- 
fect with present technology; it should however be clear that 
several other natural instances can well be conceived. Intu- 
itively speaking, it is demonstrated that due to the driving with 
thermal noise, the interaction of other modes can be effec- 
tively enhanced, giving rise to a "transistor-like" effect 0. 
We flesh out this effect at hand of two approaches following 
different approximation schemes. The first approach is es- 
sentially a weak coupling master equation, while the second 
approach makes use of stochastic samplings with respect to 
colored classical stochastic processes 1 10], which constitutes 
an interesting and practical tool to study such quantum optical 
systems of several modes in its own right. 

The observation made here adds to the insight that appears 
to be appreciated only fairly recently, in that quantum noise 
does not necessarily only give rise to heating, decoherence, 
and dissipation, providing in particular a challenge in appli- 
cations in quantum metrology and in quantum information 
science. When suitably used, quantum noise can also assist 
in processes thought to be necessarily of coherent nature, in 
noise-driven quantum phase transitions flTTl . quantum criti- 
cally fT2ll . in entanglement distillation fT3l or in quantum 
computation |[T4ll . It turns out that thermal noise, when ap- 
propriately used, can also assist in cooling. Alas, this counter- 
intuitive effect is not in contradiction to the laws of thermody- 
namics, as is plausible when viewing this set-up as a thermal 



machine or heat engine [ 15 ] operating in the quantum regime. 
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FIG. 1 . The opto-mechanical setup primarily being considered in this 
work, involving two optical modes and a mechanical one. 

The system under consideration. We consider a system of 
two optical modes at frequencies u a and cj&, respectively, that 
are coupled to a mechanical degree of freedom at frequency 
uj c . The Hamiltonian of the entire system is assumed to be 
well-approximated by H = Hq + Hi , where the free part is 
given by H = hu a a^ a+hubtf b+hu c c^ c, and the interaction 
can be cast into the form 



Hi = hg(a + 6) f (a + b)(c + c 1 "). 



(1) 



It is convenient to move to a rotating interaction picture with 
respect to Tiuj^cfia + tfb). The radiation pressure interaction 
is invariant under this transformation, while Hq simplifies to 



Hq = hAa^a - 



huJrC^C, 



where A = uj a — uj^. For most of what follows, the fre- 
quencies are chosen such that A = uo c , as we will see is 
the optimal resonance for cooling the mechanical resonator. 
This can be realized by tuning the mechanical degree of free- 
dom or the cavity mode splitting. In fact, this is exactly the 
setting proposed in Ref. |[T6l as a feasible three-mode op- 
toacoustic interaction, in an idea that can be traced back to 
studies of parametric oscillatory instability in Fabry-Perot in- 
terferometers 1 17]. Similarly, with systems of high-finesse 
optical cavities coupled to thin semi-transparent membranes 
fT8lL of double-microdisk whispering-gallery resonators fT9l 
or of opto-mechanical crystals l20l such a situation can be 
achieved. Surely numerous other architectures are well con- 
ceivable. 

In addition to this coherent dynamics, the system is as- 
sumed to undergo natural damping and decoherence - un- 
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avoidable in the opto-mechanical context. The quantum mas- 
ter equation governing the dynamics of the entire system em- 
bodying the two optical modes and the mechanical degree of 
freedom is given by 



p = Cp = p] + (£ a + £ b + C c )p, 



(2) 



with the generators being defined by C a = tzD a and 

C h = (1 + n h )nD h + n b KD b i , (3) 
C c = (1 + n c ) 7 L> c + n c7 D ct , (4) 

making use of the notation for a generator in Lindblad form 



D x (p) = 2xpx^ — {x^x,p}. 



(5) 



Here, we allow the optical bath of mode b to be in a Gibbs or 
thermal state having an arbitrary temperature. 

This type of damping reflects the plausible mechanism of 
loss. For the mechanical motion, we are primarily interested 
in the regime where u m ^> 7, such that the damping mech- 
anism of quantum Brownian motion based on some spectral 
density is virtually indistinguishable from the quantum opti- 
cal Markovian damping as for an optical mode |2H . For that 
reason, for coherence of presentation, the same type of dissi- 
pative dynamics has been chosen for the optical and mechan- 
ical modes. 

We will now discuss this given situation in two different 
pictures. The first one is a weak coupling approach leading to 
approximate analytical expressions. The second one involves 
sampling over colored classical stochastic processes. These 
methods are further discussed in the range of their validity in 
the EPAPS, where they are also compared with exact diago- 
nalization methods for small photon numbers | 23 ]. 

Description 1: Weak coupling approximation as an analyt- 
ical approach. In this approach, a picture is developed grasp- 
ing the physical situation well for small couplings g. In ad- 
dition to the actual physical baths of the three modes a, b, 
and c giving rise to dissipative dynamics, we also consider 
mode b as a further external "bath" and derive an effective 
master equation for modes a and c only. This is a good ap- 
proximation if the back action on mode b is negligible and up 
to second order in the coupling constant g. Having this pic- 
ture in mind, the Liouvillian in Eq. ^ can be decomposed as 



C = C, 



sys 



£ ba th, where ^bath = C b and 



£sys = ~^[Hq, •] + C a + £ c , £i nt = — — [Hi, •]. (6) 

Using projection operators techniques l22l . one can derive a 
master equation for the reduced system p a , c — ^b[p] 



(t-s)®p b . 



poo 

Pa A 1 ) = £sysPa,c(£)+tr 6 £ int / ds e CrS C int p a ,c 

JO 

Here C r = £ sys + £bath- Making use of the explicit expres- 
sion ([6]) for £i nt , we have 

1 f°° 

PaA 1 ) = £sysPa,c(>)-T2tr 6 [#i, / dse CrS [Hi , p a ,c (ts) £ 

fl Jo 

(J) 



In what follows, we will make a sequential approximation of 
the interaction Hamiltonian Hi and the damping mechanism. 
In order to be as transparent as possible, we mark each of the 
steps with a roman letter. 

Eq. ([7]) - up to second order expansion in the coupling g, 
which constitutes the first approximation step (a) - can also 
be written as 



1 



PaA 1 ) = £sysPa,c(£)-T^tr 6 [#i, / ds [e C - S (H t ) , p a , c (t) ®p h }} 

Ti Jo 

(8) 

where L\ acts only on the Hamiltonian Hi , corresponding to 
a "dissipative interaction picture" with respect to C r . 

We start from Eq. ([I]) and (b) neglect the term proportional 
to a^a because we assume mode a to be weakly perturbed 
from its ground state. In contrast, we allow the physical op- 
tical bath of mode b to have an arbitrary temperature and 
therefore we cannot neglect the term proportional to tfb. We 
rewrite the approximated Hi as 



H'i = Tig(a ] b + b ] a + S)(c- 



(9) 



where the operator S = tfb — n b represents the intensity fluc- 
tuations of mode b. In order to have vanishing first moments 
with respect to mode b, the mean force proportional to (b^b) 
has been subtracted, which is responsible of merely shifting 
the resonator equilibrium position. Since uj a — uj b = uj c , the 
(c) rotating wave approximation (RWA) of Eq. d9l) is 



H'( = hg(a*bc + atfJ) + hgS(c + c f ). 



(10) 



As will be explained later in more details, the first term of the 
Hamiltonian is responsible for the cooling of the mechanical 
resonator, while the second term corresponds to an additional 
heating noise. 

In order to compute the partial trace in Eq. ([5]), we need the 
two-time correlation functions of the thermal light in mode 6, 

(be c r s tf) = e~ KS n b , (Se^ s S) = e~ 2KS {nl + n b ). (11) 

The exponential functions in Eqs. §TT\ determine the time 
scale of the integral kernel in Eq. ([5]), which will be of the 
order of ft -1 . Within this time scale (d) we can neglect the 
effect of the mechanical reservoir (7 <C k), and the action of 
the map e^ s on the system operators will be 

e c ^ s c = e~^ +iuJc)s c ~ e~ iuJcS c. 

We can finally perform the integral in Eq. ([8]), and since all the 
odd moments of p b vanish, the cooling and heating terms in 
Eq. ( [TO] ) generate two independent contributions to the master 
equation, respectively 



] Pb\. 



Ik 

2 K g 2 (n 2 b +n b ) 

•Cheat — n V-^ct + L>c) , 



(12) 
(13) 
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where in calculating £heat we (e) kept only the counter- 
rotating terms. The effect of Cheat is simply a renormalization 
of the mean occupation number of the mechanical bath 



n r \-> n r 



2ng 2 (n 2 +n 6 ) 

7(4« 2 + 0;*) : 



always increasing, as expected, the effective temperature of 
the environment. Denoting with £ sys the corresponding renor- 
malized Liouvillian, the master equation can be written as 



(^sys H~ ^cool) Pa 



(14) 



With respect to Eq. ([2]), Eq. ([14]) can be numerically solved 
with much less computational resources but we have to remind 
ourselves that this approach is valid only within the RWA and 
for weak coupling: 7, # <C uo c . Another advantage of Eq. 



( 14 ) is that the corresponding adjoint equations for the number 
operators n a = a^a and n c = b^b are closed with respect to 
these operators, that is 



TLq — 2hCTIq 



((n b + l)n a - n b n c - n a n c ) , 
(n h n c - (n b + l)n a - n a n c ) + 2^n c . 



Assuming (f) that the factorization property (n a n c ) ~ 
{n a ){n c ) holds - which is essentially a mean-field approach 
which is expected to be good in case of small correlations, or, 
again as assumed, for small values of g - we can find analyti- 
cal expressions for the steady state expectation values: 



h c -r] f (n c + r/) 2 nn h n c 



2 V 4 

(h c - (n c ))7 



1/2 



where rj = 1 + n b (l + ft/7) + 2n 2 / g 2 . 

Description 2: Sampling with respect to colored station- 
ary classical stochastic processes. In this approach, we start 
from the exact dynamics Eq. ^ but treat mode b as a classi- 
cal thermal field and neglect any feed-back from the resonator. 
We substitute the bosonic operator with a complex amplitude 
b(t) ^ fit, giving rise to a semi-classical picture. The pa- 
rameter fi t can be described as a classical stochastic process 
defined by the stochastic differential equation (SDE) 



dfi t = — nfi t dt - 



~ b (dW ix) + idW iy) ), 



(15) 



with independent Wiener increments ifTOl obeying the Ito 
rules dW^dW^ = 5 a , b dt, dW^dt = 0. The dynam- 
ics of the remaining modes a and c instead, can be efficiently 
treated quantum mechanically; this is true, since for every 



single realization of the process (15), the evolution defines 
a Gaussian completely positive map and therefore the corre- 
sponding Gaussian state p a ^c\t) = £^^(Pa,c) can be de- 
scribed entirely in terms of first and second moments. The ac- 
tual quantum state of the system will in general not be exactly 



Gaussian, it can nonetheless be simulated by sampling over 
many Gaussian states associated with different realizations of 
fit'. Only the respective weight in the convex combinations 
are such that the resulting state can be non-Gaussian. The re- 
sulting state p a ,c{t) = 1&pa\c\t) will be our semi-classical 
description of the system. 

It is convenient to introduce a vector of quadratures op- 
erators u = [x c ,y c ,x a ,y a ], where xj = (j + j^)/y/2, 
Vj = ~ and j — a,c. From Eq. ([2]), we get a 

SDE for the first moments 



d{u) t 
dt 



= A t {u) t + ft 



where 
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(16) 



(17) 



can be arranged in the matrix V t 
SDE 



fit). The second moments 
= Re(uu^)t, satisfying the 



AtVt + VtAi + D + F t , 



(18) 



+ l),7(2ra c + and F t = 



dt 

where D = diag [7(271. 

ft(u)J + (u) t ff. The statistical average over many realiza- 
tion of Vt will be an estimator for the second moments of the 
quantum state V(t) = ^(V t ). In particular, the first two di- 
agonal elements give the effective phonon number of the me- 
chanical oscillator, since (n c )(t) = (Vi 5 i(£) + Vzpii) — l)/2. 
The three stochastic differential equations ( |15|16|18| ) can be 
numerically integrated in sequential order. In our simulations, 
see Fig. [2j we used the Euler method, for each time step dt 
sampling the associated Wiener increments in Eq. ( [T5] ) with 
normal distributions of variance a 2 = dt. 

Intuitive explanation of the effect of cooling by heating. 



This effect can be intuitively explained at hand of Eq. (10) 
in Description 1 : Two competing processes play here an im- 
portant role: The first term appears like a beam splitter inter- 
action between the modes a and c with a "reflectivity" given 
by the thermal fluctuations of the amplitude of mode b. This 
is responsible for the cooling of the mirror. That is to say, the 
occupation of mode b takes the role resembling the "basis of a 
transistor": A high occupation renders the interaction between 
a and c stronger, hence triggering the cooling effect. For this 
effect to be relevant, the coherence or purity of the state of b 
does not play a dominant role, and hence even thermal noise 
can give rise to cooling. This is referred to as "good noise". 
The second term corresponds to the fluctuations of the radia- 
tion pressure of mode b and it is a source of "bad" noise which 
heats the mechanical mode. 

Similarly, this effect can be studied at hand of the stochastic 



picture of Description 2, when observing Eq. ( 18 ). In addition 
to the intrinsic quantum noise described by D, stochastic fluc- 
tuations of fit generate an additional heating noise given by 
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the matrix F t . However, the same process f3 t is also contained 
in the matrix A t and corresponds to a cooling noise, up to the 
approximations identical to the above "good noise". The rea- 



son is quite evident from Eq. ( 17 ), where we observe that the 
coupling between the hot mechanical oscillator and the cold 
optical mode is mediated by the thermal fluctuations of f3 t . 
This op to-mechanical coupling, which would be zero without 
noise, leads to a sympathetic cooling of the mechanical mode. 

Example. We will now discuss the effect of cooling by heat- 
ing at hand of an example using realistic parameters in an 
opto-mechanical setting. Fig. [2] shows the effective temper- 
ature of the mechanical mode as a function of the number of 
photons in mode b: Here, effective temperature is defined as 
the temperature T of a Gibbs state 



Pc(T) 



e -huj c c^ c/{kT) 



such that (n c ) = tr(p c (T)c^c). One quite impressively en- 
counters the effect of cooling by heating, for increasing pho- 
ton number and hence effective temperature of this optical 
mode. For very large values of the photon number, the "bad 
noise" eventually becomes dominant, resulting again in a 
heating up of the mechanical mode. Note that needless to 
say, the effective temperature of the optical mode b is usually 
larger than the mechanical one by many orders of magnitude 
(approximately 10 10 K for reasonable parameters). 
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FIG. 2. Room temperature cooling with parameters reminding of 
those typical in realistic experiments [3 1: u c = 27rMHz, k = 0.2cj c , 
g = 0.3 x 10 _5 cj c , and 7 = 10~ 3 lj c . The black line shows the pre- 
dictions of the steady state using Description 1, the dots are a result 
from stochastic sampling using Description 2 (with 100 realizations), 
which qualitatively coincide well. One clearly finds that an increased 
population of mode b leads to a significant cooling of the mechani- 
cal mode - up to a point when eventually the "bad noise" becomes 
dominant. 

Summary. In this work, we have established the notion of 
cooling by heating, which means that cooling processes can be 
assisted by means of incoherent hot thermal light. We focused 
on an opto-mechanical implementation of this paradigm. We 
also introduced new theoretical tools to grasp the situation 
of driving by quantum noise, including sampling techniques 
over stochastic processes. To experimentally demonstrate this 
counterintuitive effect should be exciting in its own right. 



Putting things upside down, one could also conceive settings 
similar to the one discussed here as demonstrators of small 
heat engines |[T5l operating at the quantum mechanical level, 
where b takes the role of an "engine" and mode a of a "con- 
denser". To fully explore these implications for feasibly real- 
izing quantum thermal machines constitutes an exciting per- 
spective. It would also be interesting to fully flesh out the po- 
tential for the effect to assist in generating non-classical states 
l24l . Finally, quite intriguingly, this work may open up ways 
to think of optically cooling mechanical systems without us- 
ing lasers at all, but rather with basic, cheap LEDs emitting 
incoherent light. 
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SUPPLEMENTARY MATERIAL 

In this supplementary material, we compare the methods of 
Description 1 and 2 of the main text with an exact simulation 
of the master equation ^ in a truncated, finite-dimensional 
Hilbert space of the three involved modes, T-L = <C N (D^ 5 
<D Nc . The unique stationary state of Eq. ^ can easily be 
found numerically; a dimension d = N a NijN c of the total 
Hilbert space of, say, d < 400, is well feasible. This is obvi- 
ously an essentially exact method for small occupation num- 
bers in each of the three modes, and the error made can easily 
be estimated. This analysis, see Fig. 3, together with anal- 



ogous ones in similar regimes, shows that the methods used 
here are also suitable in the deep quantum regime. 
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FIG. 3. Stochastic simulation - introduced in Description 2 - of the 
mean number of phonons as a function of time (black line), compared 
with the exact steady state (blue line) and with the analytical approx- 
imation given in Description 1 (red line). Parameters: rib = n c — 1, 
k O.lcJc, 7 O.OlcJc, g 0.006cj c . 



